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We investigate quintessence cosmologies with a matter component consisting of particles with an 
increasing mass. While negligible in early cosmology, the appearance of a growing matter component 
has stopped the evolution of the cosmon field at a redshift around six. In turn, this has triggered 
the accelerated expansion of the Universe. We propose to associate growing matter with neutrinos. 
Then the presently observed dark energy density and its equation of state are determined by the 
neutrino mass. 



Growing observational evidence indicates a homoge- 
neous, at most slowly evolving dark energy density that 
drives an accelerated expansion of the universe since 
about six billion years [1, 2]. The origin of dark energy is 
unknown, be it a cosmological constant 01, a dynamical 
dark energy due to a scalar field (quintessence) 0,0], a 
modification of gravity [Q], or something still unexpected. 
A pressing question arises: why has the cosmological ac- 
celeration set in only in the rather recent cosmological 
past? Within quintessence models we need to explain 
a transition from the matter dominated Universe to a 
scalar field dominated Universe at a redshift z ~ 0.5. 

A similar crossover has happened earlier in the cosmo- 
logical history, namely the transition from radiation dom- 
ination to matter domination. This crossover is bound 
to happen at some time since the dilution of the energy 
density with increasing scale factor a obeys pr oc a^^ for 
radiation and pc oc for cold dark matter. At some 
moment matter must win. We suggest in this Letter that 
the presently observed crossover to a dark energy domi- 
nated Universe is of a similar type. 

We propose "Growing Matter", an unusual form of 
matter whose energy density decreases slower than the 
one of the usual cold dark matter, or even increases: 
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This may be realized by particles whose mass increases 
with time. In presence of both cold dark matter and 
growing matter a crossover to a new epoch is then nec- 
essary at some moment. In our model this transition is 
witnessed now. Similar as for the radiation-matter tran- 
sition the time for the crossover is set by the mass and 
abundance of the growing matter particles. We also spec- 
ulate that growing matter consists of neutrinos. In this 
case the abundance is computable and the crossover time 
is determined by the value of the average neutrino mass 
my. Moreover, the relation between the laboratory value 
^1/(^0 ) (for our present cosmological time to) and the 
cold dark matter density pc at the time of the crossover 
only depends on dimensionless couplings of our model. 

The appearance of a substantial growing matter com- 
ponent strongly influences the dynamical behaviour of 



the scalar field responsible for quintessence, the cosmon. 
Indeed, the possibility of a time evolution of the mass 
requires a time evolution of this cosmological scalar field. 
In our model the mass of the particles of the growing 
component obeys: 



mg{(f>) = ruge' 



f-' M 



(2) 



with M = I/vSttGtv the reduced Planck mass and rUg 
a constant. For < an increase in will induce an 
increasing mass. 

In turn, the growing matter energy density Pg influ- 
ences the evolution of the cosmon. Our approach is a 
model of "Coupled Quintessence" 0,131 • For a homoge- 
neous cosmon field the field equation [g||: 
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contains a "force" oc pg that will counteract an increase 
of (j) once f3pg is comparable to dV/d<j) ■ In our model, 
this effect will eventually dramatically slow down further 
evolution of the cosmon. For an almost static (j){t) the 
cosmon potential V{4>) will then act similar to a cosmo- 
logical constant. The expansion of the Universe therefore 
accelerates soon after 4> stops to move. The coupling be- 
tween the growing matter and the scalar field ties the 
time of onset of the accelerated expansion to the crossover 
time when Ppg becomes important. The solution of the 
"why now" problem is thus linked to the properties of 
growing matter. The mechanism we propose is similar to 
the one presented in ref. here however we suggest 

to identify the coupled matter component with the neu- 
trinos and we discuss the key role played by the growing 
matter mass. 

Let us specify our model. Besides gravity and the cos- 
mon field, for which we assume an exponential potential: 
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cosmology is determined by cold dark matter with a stan- 
dard equation of state p = 0, growing matter, baryons 
and radiation. We denote the fraction of homogeneous 
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dark energy by flh, and similarly for cold dark matter 
and growing matter by flc and fig. The cosmologically 
relevant parameters of our model are the dimensionless 
coupHngs a and /3 (eqs. I2l4p . as well as the energy den- 
sity of growing matter at some initial time, e.g. pg{teq) 
(the initial density of cold dark matter, pc{teq), can be 
translated to the present value of the Hubble parame- 
ter Hq ). We assume a flat Universe. The cosmological 
equations are the standard ones, except for the modified 
energy- momentum conservation for growing matter [^: 



3Hpg + jjPg<i> = 
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which accounts for the exchange of energy between grow- 
ing matter and the cosmon [3, 0|. In case of neutrino 
growing matter, eqs. ([3]) and |[5]) are modified by pres- 
sure terms in early cosmology. 

For the radiation and the matter dominated epochs 
in early cosmology the cosmon field follows a "tracker 
solution" or "cosmic attractor" with a constant fraction 
of early dark energy Q : 
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where n = 3 (4) for matter (radiation) . This intermediate 
attractor guarantees that the initial conditions for the 
scalar field are not fine-tuned. Observations require that 
a is large, typically a > 10 
one has 



In this "scaling regime" 
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The growing matter plays no role yet. Its relative weight 
Qg grows, however, for 7 > 1 or /? < such that growing 
matter corresponds to an unstable direction. The scaling 
regime ends once jflg has reached a value of order one. 

The future of our Universe is described by a differ- 
ent attractor 0, H] , where the scalar field and the grow- 
ing matter dominate, while baryons and cold dark mat- 
ter become negligible. The energy-momentum tensor for 
combined quintessence and growing matter is conserved 
and we define the equation of state (EOS) in the non- 
relativistic regime: 



Ph 



Ph + Pg 
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Notice that this is indeed the dark energy EOS measured 
by eg. supernovae experiments since the two coupled fiu- 
ids behave at the background level as a single conserved 
component. 



For this future attractor the expansion of the Universe 
accelerates according to (7 > 3/2): 
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For large 7 the total matter content of the Universe, 
Q.M — ilc + f^fc + f^g, will be quite small in the future, 
~ ~ 1/7- The presently observed value flm ~ 
0.25 indicates then that we are now in the middle of the 
transition from matter domination {Hm ~ 1 — S/a^) to 
a scalar field dominated cosmology {Hm ~ 1/7)- 

The limiting case 7 ^ c» admits a particularly simple 
description. In this case we encounter a sudden transition 
between the two cosmic attractors at the time tc when 
the two terms on the r.h.s. of eq. |(3]) have equal size, 
namely for aV — "Ppg or: 



(10) 



While the cosmon was evolving before this time, it sud- 
denly stops at a value (/)c = 4>{t) at tc. Thus, for t > tc 
and large 7 the cosmology is almost the same as for a 
Cosmological Constant with value V(4>c)- On the other 
hand, before tc standard CDM cosmology is only mildly 
modified by the presence of an early dark energy com- 
ponent ([6]). For large enough a this ensures compatibil- 
ity with observations of CMB anisotropies and structure 
formation. The redshift of the transition Zc may be es- 
timated by equating the potential V at the end of the 
scaling solution ([7]) to its present value. In terms of the 
present dark energy fraction Vlh.o ~ 0.75 it is given by: 
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whose solution can be approximated as 1 -I- Zc ~ 
[20/1^0^^/(3 — iVlhfi)]^^^ . In the numerical examples be- 
low we will assume a = 10 and either 7 = 5.2 or 7 = 40. 
Then we obtain numerically Zc ~ 6(5) for 7 = 5.2(40). 
Thus Zc is large enough not to affect the present super- 
novae observations. The large-7-limit is therefore com- 
patible with all present observations provided that a is 
large enough. We plot the time evolution of the different 
cosmic components and the effective equation of state for 
the combined cosmon and growing matter components in 
Fig. 1. For not too large a and 7 our model differs from 
ACDM, and we will come back below to the interesting 
possibilities of observing these deviations. 

So far we have made no assumptions about the con- 
stituents of the growing matter component. It could be 
a heavy or superheavy massive particle, say with a mass 
ITeV or lO^^TeV. Then growing matter is non-relativistic 
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at all epochs where it plays a role in cosmology. In this 
case the initial value pg(teq) has to be chosen such that 
the crossover occurs in the present cosmological epoch. 
Even more interesting, growing matter could be associ- 
ated with neutrinos. In this case our model shares cer- 
tain aspects with the "Mass Varying Neutrinos" scenario 
[l^ . although being much closer to "standard" Coupled 
Quintessence [§]. Neutrino growing matter offers the in- 
teresting perspective that no new particles (besides the 
cosmon and cold dark matter) need to be introduced. 
Furthermore, the present value of pg can be computed 
from the relic neutrino abundance and the present (aver- 
age) neutrino mass mv{to) (assuming h = 0.72): 



16eV 



(12) 



For large |/3| the neutrino mass becomes rapidly very 
small in the past such that neutrinos cannot affect the 
early structure formation. The standard cosmological 
bounds on the neutrino mass Q do not apply. 

For a given neutrino mass m^{to) our model has only 
two parameters, a and (3 (or a and 7). They will deter- 
mine the present matter density flM^to). Replacing 7 by 
^Aiito), our model has then only one more parameter, 
a, as compared to the ACDM model. For an analytical 
estimate of the relation between flMito) = l — flhitg) and 
TOi/(io) we use the observation that the ratio flg/fth (av- 
eraged) has already reached today its asymptotic value 

Mm : 
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This important relation determines the present dark en- 
ergy density by the neutrino mass and 7: 



[Phito)] 



1/4 



1.07 



V eV 
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This value will change very slowly in the future since the 
value 7 = 5.22(800) for the maximal (minimal) neutrino 
mass (no sterile neutrinos) m,y{to) = 2.3eT^(0.015ey) 
must indeed be large and w is therefore close to —1, cf. 
Fig. 1. The late dark energy density is essentially deter- 
mined as the neutrino energy density times 7. Its actual 
value is given by the value of the scalar potential at the 
crossing time ic, i-e. 9M^H'^{tc)/2a'^. Since the equation 
of state im is today already near the asymptotic value 
([9]), cf. Fig. 1, we can relate it to the neutrino mass 
(nh.,0 ~ 3/4) by eqs. dgil 



-1 + 
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This remarkable expression yields m^ito) < 2.4eV for 
w < -0.8. 



How can our model based on growing matter be tested 
and constrained? First of all, the presence of early 
dark energy manifests itself by the detailed peak location 
of the CMB anisotropies |13|. the change in the Hnear 
growth of cosmic structures |l4l.[lB|. and the abundance 
and properties of nonlinear structures [l^ . 

Second, for not too large 7 there would be a sizable 
fraction of growing matter today (for neutrino growing 
matter this would require rather large neutrino masses). 
Then the present matter density pM — Pc + Pb + Pg differs 
from the (rescaled) matter density in the early Universe 
Pc + Pb- This may affect the matching of the present val- 
ues of flM and flb/flM obtained from supernovae, baryon 
acoustic oscillations and clusters, with determinations 
from the CMB at high redshift, through the value of teq 
and the baryon content of the Universe at last scattering. 
This effect is small for large values of 7 (small neutrino 
mass). 

Third, growing matter can affect the formation of 
structures in the late stages. For very massive parti- 
cles, growing matter would consist of relatively few par- 
ticles which have presumably fallen into the cold dark 
matter structures formed in early cosmology. For scales 
smaller than the range of the cosmon interaction these 
particles feel a strong mutual attraction, enhanced by a 
factor (2/?^ -I- 1) as compared to gravity. This force is 
mediated by the cosmon 0, [l3|- Thus, once a sufficient 
fig is reached, the growing matter structures Spg grow 
rapidly. They will influence, in turn, the structures in 
baryons and cold dark matter once the gravitational po- 
tential of the growing matter structures becomes compa- 
rable to the one of the cold dark matter structures. This 
happens rather late, especially for large 7 since growing 
matter constitutes only a small fraction of the present 
matter density in this case. 

The condition for the onset of an enhanced growth of 
Spg requires that the average cosmon force ^ 20^^g is 
comparable to the average gravitational force ~ ^m- 
This happens flrst at a redshift z^g somewhat larger 
that the crossover Zc (eq. [TT]). At this time the 
scaling solution is still vaHd, with f^M ~ \ — "ijo? 
and ^g{z) = [(l + zc)/{l + z)f^-^^ng{zc), ng{zc) « 
-fflvizc) ~ 37/(20;^), resulting in: 
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For large 7 one finds Zeg quite close to Zc such that the en- 
hanced growth concerns only the very last growth epoch 
before the accelerated expansion reduces further Hnear 
growth in the dark matter component. For heavy grow- 
ing matter this results in an enhancement of as as com- 
pared to the ACDM model, which may be compensated 
by a slower growth rate before Zc due to early dark energy 

3, El- 

We notice that the fluctuation growth rate / = 
of the dominant form of matter can be estimated ana- 
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lytically during both the scaling phase and the future 
attractor. During the first phase the matter fiuctuations 



grow as 



h = -(-1 + V240m + 1) 



(17) 



and are therefore slowed down with respect to the stan- 
dard matter dominated rate; for large a the modification 
is small, /i « 1 — 9a^^/5. When the growing component 
becomes dominant, its growth rate can be evaluated as 
a function of qv7 that for a ^ \ and 7^1 reduces to 
/2 « 1.04a^/7 [1, [iBl- Since in our model we are not yet 
on the final attractor the present growth rate of the dark 
matter component will be quite smaller than /2. It is 
interesting to remark that in principle an estimate of /i 
and /2 would completely fix both parameters a, 7. 

An enhanced growth of 5pg concerns only structures 
with size smaller than the range of the cosmon- 
mediated interaction. In a cosmological situation we have 
to solve, in principle, the coupled system of linear fiuctu- 
ation equations in 0, pg, pc, etc... For a rough estimate 
of the (density dependent) cosmon mass we consider a 
fixed pg. 
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During the scaling solution |(7]) the cosmon range is given 
by 0: 



(19) 



After tc the evolution of 4> essentially stops, resulting in 
a constant range = ^/2/{3H{tc)). 

A different regime of growth appHes for I > l^. A 
"window of adiabatic fiuctuations" opens up in the range 
Iff, < I < H^^ where the fiuctuations of the coupled cos- 
mon fiuid and growing matter can be approximated as a 
single fiuid. In this regime the enhanced growth is weak- 
ened by the small range of the cosmon interaction. 

Neutrino growing matter was relativistic in earlier 
time, so that free streaming prevents clustering. For 
/? < neutrinos have actually remained relativistic much 
longer than neutrinos with constant mass. In the limit 
of large 7 one can estimate that the neutrinos are rela- 
tivistic at a < an where 



aR « {mu{tQ)/Tufl) 



-1/4 



0.11 



leV 



-1/4 



(20) 



which corresponds to £ (2— 10) for m^^o G (0.015—2.3) 
(we use mjj[aii) = T^io-n), and approximate pg ^ const, 
so that nil, ~ Va^, V ^ const, Ti, = T^fi/a where T^^ is 
the present neutrino temperature for massless neutrinos) . 
The growth of neutrino structures only starts for z < 
zr. Even then, neutrinos cannot cluster on scales smaller 
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Figure 1: Cosmological evolution for neutrino growing matter 
for Q = 10,/3 = —52 and m^fi = 2.3eV. Panel a): density 
fractions firad (black, dots), fie (dark green, dot-dashes), fib 
(blue, short dashes), ilh (red, long dashes) , ilg (light green, 
solid). Panel b): blow out of panel a) near the present time. 
Panel c): total equation of state w^ff = ptot/ptot (red, long 
dashes); combined EOS of cosmon and neutrinos (blue, short 
dashes); and EOS of cosmon alone (green, solid). Panels a) — 

c) remains almost identical for heavy growing matter. Panel 

d) : neutrino energy density (red, long dashes), neutrino mass 
(green, solid) normalized to unity today. The dotted curve 
represents the energy density of always non-relativistic heavy 
growing matter. 
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than their "free streaming scale" Ifs . This scale is given 
by the time when the neutrinos become non-relativistic, 
eq. {201), close to H-^ur) « 200 (m^/leF)^/^ ft-iMpc 
e (100- 1000) h-^Mpc. 

On the other hand, the range of the effective cos- 
mon interaction increases, eq. (fTOj) until it reaches l^. 
For scales within the window Ifs < I < l^p the neu- 
trino clustering is strongly enhanced (for z < z^g) due 
to the additional attractive force mediated by cosmon 
exchange. This enhanced clustering starts first for scales 
close io Ifs ■ One may thus investigate the possible for- 
mation of lumps with a characteristic scale around Ifs- 
For the range < I < one expects again an adia- 
batic growth of the coupled neutrinos and cosmon fiuc- 
tuations, approximated by a single fiuid. In summary, on 
large scales I > Ifs the neutrino fiuctuations grow similar 
to the heavy growing matter fluctuations. The growth 
starts, however, only very late for z > zr and only from 
a low level given by the tiny fluctuations in a relativistic 
fluid at Zr. Furthermore, neutrino fluctuations with a 
scale I < Ifs are suppressed by free streaming. 

In this context one may ask to what extent the local 
variations of the cosmon fleld (fi affect the local values of 
the neutrino mass and therefore the relation of cosmol- 
ogy to the possible outcome of laboratory experiments. 
Such variations in (j) are induced by local neutrino con- 
centrations. For neutrino lumps the perturbations of (j) 
are of the order of /3 times the Newtonian potential $. 
Thus one expects local variations of the growing matter 
mass of the order of Amg/nig « In all observable 

(non compact) astrophysical objects one has $ < lO^'' 
and therefore it turns out that for /? < 100 the spatial 
variations of nig should not be very large. In particu- 
lar, we can safely identify the cosmological value of the 
neutrino mass m^{to) with the one measured on Earth. 

A close look at Fig. 1 shows oscillations of fig, start- 
ing around Zc and being damped subsequently. Both the 
oscillation period and the damping time can be under- 
stood in terms of the eigenvalues of the stability matrix 
for small fluctuations around the future attractor solu- 
tion 0, We note, however, that the oscillations con- 
cern only the relative distribution between fig and flh, 
while the sum flh + ^g = I ~ i^M remains quite smooth. 
A detection of the oscillations by investigations of the 
background evolution, like supernovae, seems extremely 
hard - the luminosity distance is a very smooth func- 
tion of z. For neutrino growing matter the oscillations of 
m,^{t) around its average value induce an uncertainty in 
the estimate of the relation between 51g(io) and m,y{to) 
(fT3l) ranging from 25% for 7 = 5.2 (m^ito) = 2.3eV) to 
100% for 7 = 40 (m^(to) = 0.3eV ). The relations (fTMS]) 
all involve the averaged neutrino mass. 

While for large enough 7 and a the cosmology seems 
rather reaHstic, one may ask if our proposal for a res- 
olution of the "why now" problem has not introduced 
other unnaturally small parameters. Indeed, extrapolat- 



ing the masses of the growing matter particles back to the 
Planck time may result in extremely small masses if 
is large, given typical values (j){tpi) ~ 0, a(t>{tQ)/M « 276 
(from eq. [71). However, we have explored here only the 
simplest possibility of constant a and (3. It is well con- 
ceivable (and quite likely) that in a fundamental theory 
a and (3 are functions of 4>. Slow changes will not affect 
our phenomenological discussion which only concerns a 
rather small range of (f)/M. In contrast, extrapolations 
back to the Planck epoch or Inflation could look com- 
pletely different. Our scenario does not need a huge over- 
all change of the mass of the growing matter particles. 
For neutrinos a growth of the mass by a factor 10'', cor- 
responding in the seesaw mechanism to a decrease of the 
right handed neutrino mass from M to lO^^GeT^ , would 
largely be sufficient, provided a fast change happens dur- 
ing recent cosmology. We also have assumed here that 
cold dark matter has a negligible coupling to the cosmon. 
While the cosmon coupHng to baryons must be very small 
in order to remain compatible with the tests of General 
Relativity, the coupling between cold dark matter and 
the cosmon is only restricted by cosmology and needs 
exploration in our context. 

The most crucial observational issues can be under- 
stood by concentrating on constant parameters a, (3 (and 
possibly a constant cosmon-cold dark matter coupling). 
It will be a challenge to measure them or to falsify the 
growing matter scenario. For neutrino growing matter a 
determination of a and (3 would fix the neutrino mass, 
allowing for an independent test of this hypothesis by 
comparing with laboratory experiments. 



[1] P. Astier et a l, Astroii. Astroph ys. 447 31 (2006); T. M. 

Davis et al, astro-ph/0701510 
[2] D.N. Spergel et al (WMAP collaboration) ApJS, 170, 

377 (2007) arXiv:astro-pli/0603449 
[3] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 
[4] C. Wetterich, Nucl. Phys. B 302 (1988) 668. 
[5] P. J. E. Peebles and B. Ratra, Astrophys. J. 325 (1988) 

L17. 

[6] G.R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 

484, 112 (2000) 
[7] C. Wetterich, Astron. Astrophys. 301 (1995) 321 

|arXiv:hep-th/94 08025|. 
[8] L. Amendola, Phys. Rev. D62 (2000) 043511 

fastro-ph/9908023 1, L. Amendola, D. Tocchini-Valentini, 

Phys. Rev. D66 (2002) 043528 |astro-ph/0111535 1. 
[9] C. Wetterich, Nucl. Phys. B 302 (1988) 645. 
[10] G. Huey & B. Wandelt, Phys. Rev. D74 (2006) 023519 
[11] M. Doran, G. Robbers and C. Wetterich, Phys. Rev. D 

75 (2007) 023003 |arXiv:astro-ph/0609814 1. 
[12] P. Gu, X. Wang and X. Zhang, Phys. Rev. D68, 087301 

(2003); R. Pardon, A.E. Nelson and N. Weiner, JCAP 

0410 (2004) 005; A. W. Brookfield et al., Phys. Rev. 

Lett. 96 (2006) 061301; N. Afshordi, M. Zaldarriaga & K. 

Kohri Phys. Rev. D72 (2005) 065024; O. E. Bjaelde et al. 



6 



arXiv:0705.2018; K. Ichiki & Y. Keum, arXiv:0705.2134 
[13] M. Doran, M. J. Lilley, J. Schwin dt and C. W etterich, 

Astrophys. J. 559 (2001) 501 |arXi v:astro-ph/00 12139 1. 
[14] P. G. Ferreira and M. Joyce, Phys. Rev. Lett. 79 (1997) 

4740 |arXiv:astro-ph/9707286|. 
[15] M. Doran, J. M. Schwindt and C. Wetterich, Phys. Rev. 



D 64 (2001) 123520 |arXiv:astro-ph/0107525 1. 
[16] M. Bartelmann, M. Doran and C. Wetterich, Astron. & 

Astrophys. 454 (2006) 27 |arXiv:astro-ph/0507257|. 
[17] C. Wetterich JCAP 0310 (2003) 002 [arXiv: 

|hep-ph/0203266] 



